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The pointing accuracy needed to support Magellan ’s Synthetic Aperture Radar 
mapping of Venus places stringent requirements on navigation accuracy. This need 
is met with a combination of two-way Doppler and narrowband A Very Long Base- 
line Interferometer (AVLBI) data, which are capable of determining the spacecraft’s 
orbit to the required level, typically about one-kilometer position uncertainty. Dif- 
ferenced Doppler (two-way Doppler minus three-way Doppler) is also capable of 
meeting mission navigation requirements, and serves as a backup to narrowband 
i, - AVLBI. The Magellan Project specifies that the turn-around time for processing 

* narrowband AVLBI data must be 12 hours or less, a very difficult requirement 

to meet operationally. In this article, the use of phase-delay data, taken from a 
Connected-Element Interferometer (CEI) with a 21-km baseline, for Magellan orbit 
determination was investigated to determine if navigation performance comparable 
with narrowband AVLBI and differenced Doppler could be achieved. CEI possesses 
an operational advantage over AVLBI data in that the observables are constructed 
in near-real time, thus greatly reducing the turn-around time needed to process the 
data, relative to the off-line system used to generate AVLBI observables. Unfortu- 
nately, the results indicate that CEI data are much less powerful than narrowband 
AVLBI and differenced Doppler for orbiter navigation, although there was some 
marginal improvement over the navigation performance obtained when only two- 
way Doppler data were used. 


I. Introduction 

Magellan will arrive at Venus in August 1990, and 
after a short period of time devoted to orbit trim and sys- 
tems checkout, the spacecraft will begin its 243-day pri- 
mary mission phase, the mapping of Venus. Navigation of 
Magellan will be very challenging due to the high pointing 
accuracy required by the Synthetic Aperture Radar (SAR) 
system. The SAR pointing accuracy requirements corre- 
spond to a navigation accuracy of about a one-kilometer 


(lcr) position uncertainty at periapsis. Furthermore, this 
level of performance must be maintained over most of the 
primary mission, except for an 18-day period around su- 
perior conjunction, during which mapping operations are 
suspended. 1 


1 S. Mohan, et al., “Magellan Navigation Plan,” JPL D-14S0 (inter- 
nal document), Revision B, Jet Propulsion Laboratory, Pasadena, 
California, March 23, 1988. 
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Navigation studies carried out early in Magellan’s 
history 2 discovered that Doppler-based orbit determina- 
tion, which had been sufficient for previous missions such 
as Mariner-Mars ’71 and Viking, would not meet naviga- 
tion requirements and that new radio metrics, such as nar- 
rowband A Very Long Baseline Interferometry (AVLBI) 
and differenced Doppler, would be needed [1]. It is now 
well known that the orbit determination accuracy obtained 
with two-way Doppler is dependent on the orbit orienta- 
tion with respect to the plane of the sky, which is the 
plane perpendicular to the Earth-spacecraft line of sight 
[2,3]. Recent studies done in support of Magellan navi- 
gation planning 1 [4] show that the combination of either 
narrowband AVLBI or differenced Doppler with two-way 
Doppler produces improvements of a factor of 2 to 10 in 
orbit determination uncertainty over Doppler-only naviga- 
tion. 

The major drawback involved with the use of nar- 
rowband AVLBI data is the time required to process the 
raw measurement data into observables. The Magellan 
Project requires a turn-around time of 12 hours for each of 
the two narrowband AVLBI scans scheduled daily during 
the primary mission. Connected Element Interferometry 
(CEI), a new technique currently under development, can 
be used to generate medium accuracy (100-400 nrad for in- 
dividual measurements) angular observations in near-real 
time, using relatively close (10-100 km apart) stations [5]. 
This study investigated the possibility that CEI angular 
measurements could be used to infer angular rate informa- 
tion, which is directly observed with narrowband AVLBI. 
It was hoped that this approach might be capable of pro- 
ducing orbit determination accuracy comparable with that 
of AVLBI, thus offering an alternative data type for Mag- 
ellan navigation that would reduce the operational strain 
imposed on the Deep Space Network (DSN). 

II. The Magellan Venus Mapping Orbit 

Magellan orbital elements relative to Venus’ equator 
are given in Table 1. Since Venus has a nearly spherical 
gravity field, the spacecraft’s orbit remains nearly fixed 
in orientation with respect to Venus throughout the pri- 
mary mission. The mapping orbit’s orientation with re- 
spect to the plane of the sky does, however, change during 
the primary mission due to the relative motion of Venus 
and the Earth. The orbit geometry chosen for analysis 
in this study represents Magellan’s orbit as seen from the 
Earth on December 4, 1990, 120 days after Venus Orbit 


2 J. Ellis and R. Russell, “Earth-Based Determination of a Near- 
Circular Orbit About a Distant Planet,” JFL Technical Memo- 

randum 391-406 (internal document), Jet Propulsion Laboratory, 
Pasadena, California, March 9, 1973. 


Insertion (VOI+120). Table 2 gives the orbital elements 
in the plane-of-the-sky coordinate system for this date. 

HI. Tracking Schedule 

Nominally, trajectory solutions for Magellan are gen- 
erated daily, based on 24 hours of continuous two-way 
Doppler data from Deep Space Stations (DSSs) 14, 42, and 
61 and on two 10- minute AVLBI scans using the baselines 
defined by DSS 14 to 42 and DSS 14 to 61, taking place 
during the same 24-hour period. Doppler data are not 
taken for a 1-hour period centered about each periapsis, 
since the spacecraft’s antenna is pointed at Venus during 
that time. Doppler data are also not taken from 10 min- 
utes before to 5 minutes after each apoapsis while a star 
sighting is being done to update the spacecraft’s inertial 
attitude reference system. 3 Since Magellan’s orbit period 
is 3.15 hours, the data going into each trajectory solution 
are spread over just under eight revolutions. There are cer- 
tain time periods within the primary mission in which the 
declination of the spacecraft is too low for it to be viewed 
from the DSS 14 to DSS 61 baseline, thus leaving out one 
of the two daily AVLBI scans. The epoch for this investi- 
gation, December 4, 1990, was intentionally chosen during 
one of these periods in which only one AVLBI baseline 
is available, since it is these periods in which the mission 
navigation requirements can only marginally be met. 

Simulated CEI data were generated using the DSS 13 
to DSS 15 baseline within the Goldstone DSN complex. 
This baseline is 21 kilometers in length and is oriented in 
a near north-south direction. During the 24-hour period 
used in the study, the spacecraft was in view of DSS 13 
to DSS 15 for just under three revolutions, with a total 
view time of six hours available, subject to the restrictions 
described above regarding time periods during which the 
spacecraft could not be tracked. Two different CEI data 
rates were used in the CEI navigation simulations, either 
5 points/hour or 10 points/hour. 

One case was studied which included differenced 
Doppler data (two-way Doppler minus three-way Doppler) 
in place of narrowband AVLBI data. The data schedule 
for this case contained the Doppler data from the baseline 
case plus one pass (approximately two hours) of differenced 
Doppler data using the DSS 14 to DSS 42 baseline. Dif- 
ferenced Doppler data are constructed by differencing the 
Doppler signal received at a Deep Space Station operating 
in a coherent (two-way) mode with the Doppler signal re- 
ceived at another, distant station (known as "three-way” 


3 S. Mohan, et ah, op. cit. 


Doppler). In theory, the information content of difFerenced 
Doppler data is equivalent to that of narrowband AVLBI 
data. In practice, the narrowband AVLBI data are more 
accurate due to the cancellation of the effects of several 
error sources common to both receiving stations through 
the use of a quasar (extragalactic radio source, or EGRS). 
Because of this, it is necessary to obtain more differenced 
Doppler data than narrowband AVLBI data to get com- 
parable navigation performance. 


IV. Filter Model 


D. Data Accuracy 

(1) Two-way Doppler: 10.0 mm/sec (60-sec count 

time). 

(2) Narrowband VLBI: 0.13 mm/sec (equivalent to 
1 mHz). 

(3) Differenced Doppler: 0.39 mm/sec (equivalent to 
3 mHz). 

(4) CEI phase-delay: 1 mm (50 nrad, 4 mm (200 
nrad), or 8 mm (400 nrad). 


A. Type 

The filter used was a batch-weighted, least-squares, 
with the exception that narrowband AVLBI and differ- 
enced Doppler data were processed sequentially due to the 
presence of stochastic parameters, as noted below. 

B. Estimated Parameters 

(1) Spacecraft epoch state vector. 

(2) Venus harmonic gravity coefficients up to degree 
3, order 3. 

(3) Spacecraft drag coefficient. 

(4) Stochastic station oscillator frequency error— this 
parameter is only included for cases with narrow- 
band AVLBI data. 

(5) Stochastic troposphere error— this parameter is 
also included only for cases with narrowband 
AVLBI data. 

(6) Stochastic Doppler frequency error— this param- 
eter is only included for the cases with differenced 
Doppler data. 

C. Considered Parameters 

(1) Venus harmonic gravity coefficients of degree and 
order 4-5. 

(2) Mass of Venus. 

(3) Seven mass concentrations ( mascons ), see Table 
3 for a list detailing the name, location, size, and 
a priori uncertainty of each. 

(4) Station coordinates of DSSs 14 (Goldstone), 42 
(Canberra), and 61 (Madrid). 

(5) Impulsive maneuver representing daily momen- 
tum wheel desaturation. 

(6) Quasar location uncertainty (right ascension and 
declination). 


V. A Priori 1 Uncertainties 

A. Estimated Parameters 

(1) Spacecraft epoch state: 10 6 km in position and 
1 km/sec in velocity. 

(2) Venus 3x3 gravity harmonics. 4 

(3) Spacecraft drag coefficient: 0.5 (50-percent un- 
certainty). 

(4) Stochastic station oscillator frequency error: the 
t t of this parameter is zero except during the 10- 
min AVLBI scan, when the a = 0.3538 mHz, 
with a time constant r = 27.6 sec. 

(5) Stochastic troposphere error: this parameter also 
has a <r of zero except during the AVLBI scan, 
during which the a = 1 223 mHz, with a time 
constant r = 60.7 sec. 

(6) Stochastic Doppler frequency error: the a is zero 
except during the differenced Doppler pass, dur- 
ing which it has a <r — 0.55 mHz, with a time 
constant r = 3604 sec. 


B. Considered Parameters 

(1) Remaining Venus gravity harmonics up to 5 x 5 4 

(2) GM of Venus: 1.0 km 3 / sec 2 . 


(3) Mascons: see Table 3. 

(4) Station coordinates for DSSs 14, 42, and 61: 


Coordinate 
Spin rad., m 
Z-height, m 
Long., deg 


14 

39 

54 

0.74 

0.81 

0.79 

8.70 

8.70 

8.70 

1.8 x 10~ 5 

2.1 x 10“ 5 

2.2 x 10 


(5) Momentum wheel desaturation burn: 8.3 x 10 
km/sec. 


4 Ibid., Appendix E. 
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(6) Quasar location uncertainty: 167 nrad in both 
right ascension and declination. 

VI. Results 

A total of eight simulation cases was run, all using the 
Magellan delivery of the Orbit Determination Program on 
a Sun 3/260 workstation. Table 4 contains a summary 
of the results, with navigation performance being defined 
in terms of the lcr uncertainty in the longitude of the as- 
cending node in the plane of the sky, the dominant orbit- 
determination error. Figure 1 shows a graphic comparison 
of the eight cases. This figure shows that, of the five cases 
run using simulated CEI data, only one yielded results 
that were appreciably better than the Doppler-only case. 
Recalling that only one scan of narrowband AVLBI data 
was included in the AVLBI case (taken from the DSS 14 
to DSS 42 baseline; the spacecraft cannot be viewed from 
the DSS 14 to DSS 61 baseline at this time), Fig. 1 shows 
that the performance obtained using two hours of differ- 
enced Doppler data was actually better than that obtained 
with AVLBI data. This is not always the case however, 
when narrowband AVLBI scans can be taken using both 
the DSS 14 to DSS 42 and the DSS 14 to DSS 61 baselines. 

A breakdown of the root sum square (RSS) contribu- 
tion of different error source groups to the total RSS error 


in longitude of the ascending node is given in Table 5. Fig- 
ure 2 shows a graphic version of the data in Table 5, the 
four cases in the figure corresponding to the four cases in 
the table. Figure 2 shows clearly that the CEI data are 
much more sensitive to gravity field uncertainty and mas- 
cons than the AVLBI data. Even excluding all consider 
error sources, the comparison of the computed (state-only) 
solution uncertainty shows that narrowband AVLBI data 
are much more powerful than CEI data in this situation. 
Further, Fig. 1 shows that, even if the CEI data were ac- 
curate to 1 millimeter (which is very optimistic), the sen- 
sitivity of this data type to the consider errors becomes 
extreme. 

VII. Conclusions 

An error analysis was conducted to determine the use- 
fulness of CEI phase-delay data in the navigation of the 
Magellan mission. The results indicate that CEI data are 
much less capable than narrowband AVLBI and differ- 
enced Doppler data in this situation, both in information 
content, which is reflected in the comparison of computed 
solution uncertainties, and in sensitivity to consider error 
sources such as gravity field uncertainty. Based on the re- 
sults presented in this article, this data type would likely 
not be capable of meeting Magellan navigation require- 
ments during any portion of the primary Venus mapping 
mission. 
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Table 1. Magellan orbital elements relative to Venus equator 


Epoch: December 4, 1990 18:20:42 UTC 
Semi-major axis: 

a 

10189.7 

Eccentricity: 

e 

0.38168 

Inclination: 

* 

85.9743° 

Argument of peri apsis: 

w 

170.016° 

Longitude of ascending node: 

U 

354.205° 

Altitude at periapsis: 


249.50 km 

Altitude at apoapsis: 


8027.0 km 

Sun- Earth-probe angle: 


<£> 

o 

o 


Table 2. Magellan orbital elements relative to the plane of the sky 


Epoch: December 5, 1990 22:42:30 UTC 


Semi-major axis: 

a 

10189.7 km 

Eccentricity: 

e 

0.38168 

Inclination: 

i 

55.7623° 

Argument of periapsis: 

w 

84.2673° 

Longitude of ascending node: 

Q 

10.0999° 


Table 3. Mascons considered in Magellan orbit determination 1 


Name 

Lat, deg 

Long, deg 

Radius, km 

a km 3 /sec 2 

Ovda Regio 

— 5.4 

95.6 

900 

0.040 

Leda 

42.5 

55.2 

600 

0.012 

Hestia 

6.4 

64.6 

350 

0.010 

Tellus 

34.0 

83.2 

900 

0.020 

Bell Regio 

28.9 

49.0 

500 

0.010 

FN27 

7.1 

78.0 

350 

0.010 

Maxwell 

65.0 

10.0 

900 

0.040 

» All seven mascons lxave a radial position of 6000 km; for comparison, Venus has 
a mean radius of 6052 km. 


Table 4. Summary of results for all cases 


Case 

Quasar-S/C Angle, deg 

Node Uncertainty, 1 a 

CEI 50 nrad, 10 pt/hr 

26 

0.25650° 

CEI 200 nrad, 10 pt/hr 

26 

0.10530° 

CEI 400 nrad, 10 pt/hr 

26 

0.10140° 

CEI 200 nrad, 5 pt/hr 

26 

0.09959° 

CEI 200 mad, 10 pt/lir 

3 

0.08809° 

Doppler only 

- 

0.10450° 

Doppler -f AVLBI 

26 

0.03679° 

Doppler + differenced Doppler 

- 

0.02963° 


Table 5. Root-sum-square error breakdown of node error for selected cases 



RSS Contribution, deg a 


Error Source 

Case 1 

Case 2 

Case 3 

Case 4 

Computed 

0.0242 

0.0479 

0.0500 

0.0532 

Gravity field 

0.0195 

0.0650 

0.0700 

0.0801 

Mascons 

0.0071 

0.0229 

0.0246 

0.0281 

Quasar location 

0.0183 

0.0265 

0.0555 

0.0158 

Station locations 

0.0003 

0.0001 

0.0001 

0.0001 

Momentum wheel desat 

0.0002 

0.0008 

0.0009 

0.0010 

RSS total, deg 

0.0368 

0.0880 

0.1053 

0.1014 

a Case 1: Baseline operational case with Doppler, narrowband AVLBI 


Case 2: Doppler, 4-mm 

CEI, data rate 10 pt/hr, quasar-S/C angle 3 


Case 3: Doppler, 4-mm 

CEI, data rate 10 pt/hr, quasar 

-S/C angle 26 

0 

Case 4: Doppler, 8-mm 

CEI, data rate 10 pt/hr, quasar^S/C angle 26 

0 


1-mm CEi, 10/hr 


4-mm CEi, 10/hr 
8-mm CEI, 10/hr 
4-mm CEI, 5/hr 
4-mm CEI, 3 deg sep 
DOPPLER ONLY 
DOPPLER + DDOP 
DOPPLER + VLBl 


Fig. 1. Navigation performance obtained with different data strategies. 




0.04 


0.08 


0.12 0.16 
StGMA NODE, deg 


0.20 


0.24 


0.28 



Fig. 2. RSS error breakdown for selected data strategies. 



